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Abstract 

The deconfinement phase transition is studied in the ensemble canonical 
with respect to triality. Since this ensemble implies a projection to the zero 
triality sector of the theory we introduce a quantity which is insensitive to 
Z(N C ) symmetry but can reveal a critical behaviour in the theory with dy- 
namical quarks. Further, we argue that in the canonical ensemble description 
of full QCD there exist domains of different Z(N C ) phases which are degener- 
ate and possess normal physical properties. This contradicts the predictions of 
the grand canonical ensemble. We propose a new order parameter to test the 
realization of the discrete Z(N C ) symmetry at finite temperature and calculate 
it for the case of Z(2) gauge fields coupled to fundamental fermions. 



Work supported by Bundesministerium fiir Wissenschaft, Forschung und Kunst of Austria. 



1 Introduction 



Ever since the end of the seventies there has been increasing interest towards dis- 
crete global and local gauge symmetries. In 1978 G. 't Hooft fT] and G. Mack || 
conjectured that the discrete center Z(N C ) of an underlying gauge group G can be of 
crucial importance for quark confinement. Two years later the Z(N C ) mechanism of 
confinement by means of vortex condensation was presented ||. The next important 
step in understanding the role of discrete symmetries was the awareness of the deep 
connection between the spontaneous breaking of the global Z(N C ) symmetry and 
the deconfinement phase transition in pure gauge models j|. Recently, different 
aspects of global and local discrete symmetries have been an object of numerous 
investigations both on the lattice and in continuum. This article contributes to pre- 
vious studies from the point of view of the canonical ensemble (CE) description of 
hot gauge theories with fundamental quarks. 

We propose a new order parameter to test a realization of the discrete Z(N C ) 
symmetry at finite temperature and calculate it for the case of Z(2) gauge fields 
coupled to fundamental fermions. It is argued that domains of different N-ality in 
full QCD may exist, are degenerate and possess normal physical properties. What 
we would like to stress here is that in some important aspects the two possible 
descriptions produce different results and, consequently it would be safer to use the 
CE with respect to triality. 

We begin over viewing the standard Z(N C ) picture in QCD at finite temperature 
[|]]. Gluon fields are strictly periodic in time, while quarks are anti periodic with 
a period being the inverse temperature. Pure gauge theory has an exact Z(N C ) 
global symmetry. The gauge invariant operator, the Polyakov loop (PL), Lg = 
jfTr II^i Uo(x,t) transforms under Z(N C ) global transformations as 

2ni 

Lg -> ZLg, Z = exp[— n], n = 0, N c - 1. (1) 

The PL can be used as an order parameter to test Z(N C ) symmetry in pure gauge 
theory. The expectation value of the PL is interpreted as the free energy of a probe 
quark F q immersed in a pure gluonic bath 

<L*>=exp(~F ff ). (2) 

Unbroken Z(N C ) symmetry implies < L >= and F q = oo. When the global 
Z(N C ) symmetry is spontaneously broken, < L > differs from zero and F q has a 
finite value, i.e. it costs only finite energy to create a single quark in the gluonic 
bath. To our knowledge, the first objection against this lore was discussed in || 
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where it was pointed out that in a system with a finite UV cutoff the free energy 
should not diverge. Usually, MC and analytical calculations are performed also with 
periodic boundary conditions in space directions. It has been shown, however, that 
Gauss' law and this p.b.c. in space are inconsistent unless the sum of quark and 
gluon colour charges vanishes, which is impossible since they have different triality. 
It was concluded that for space p.b.c. the expectation value of the PL is not the 
free energy of a heavy quark. The second trouble appears if we realize that in the 
spontaneously broken phase < L > may pick up N c different values corresponding to 
N c equivalent minima of the free energy. Thus, L can be negative or even complex. 
Eq.(fJ) tells us that the free energy could be a complex number. This gives rise to 
doubts that Eq.@ has the proper physical interpretation. When dynamical quarks 
are included the picture becomes more complicated and new troubles appear. The 
fermion determinant generates loops going around the lattice a number of times 
which is not a multiple of N c . Such loops present a propagation of single quarks 
and transform non trivially under Z(N C ). This means that dynamical quarks break 
Z(N C ) symmetry explicitly and screen sources of heavy quarks at any temperature. 
The expectation value of the PL prefers the phase with argL = which provides the 
minimum of the free energy. Other Z(N C ) phases with argL = ^im, k — 1, . . . , N c — 1 
become metastable. They possess, however, such unphysical properties as a complex 
free energy or entropy ||. Recently, it has been discovered that chiral symmetry 
is not restored in Z(N C ) phases 0. This led to a reexamination of degenerate 
Z(N C ) phases and interfaces between them in pure gauge theory || §. Although 
not common but the most popular opinion is that all Z(N C ) phases correspond to 
the same physical state and the interface is unphysical. 

These properties reflect the GCE description of QCD with respect to triality. In 
our previous papers we introduced an ensemble canonical with respect to triality [[10 



(see also pl| j and studied the QCD phase structure within this ensemble |L2], [T 
This ensemble reveals the following important properties of QCD with dynamical 
quarks: 

1) In the low temperature phase every state has triality zero. In the deconfine- 
ment phase the system as a whole has zero triality. Since all Z(N C ) noninvariant 
variables are projected out of the theory, the PL itself has only little meaning in 
the CE. A single quark does not appear in the spectrum. Since the Debye mass is 
nonzero, chromoelectric fields of quarks are screened. To distinguish the low temper- 
ature phase from a non-confining phase one can use the phase of the SU(N C )/Z(N C ) 
PL. Since this phase is insensitive to Z(N C ) transformations it can be used also in 
the theory with dynamical quarks ]Tj| [H| . The deconfinement phase transition can 
be seen in the CE similarly to pure gauge theory |13 . 

2) Metastable minima with unphysical properties are absent, all Z(N C ) phases are 
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degenerate just as in a pure gauge model. It follows, Z(N C ) phases and corresponding 
domain walls may exist in full QCD JlO, 0, 13 . 



3) What is not so obvious is the behaviour of the quark condensate in different 
Z(N C ) phases. It was shown, for example, that chiral symmetry is always broken 
in the metastable minimum of SU(2) theory and the quark condensate differs from 
zero in this minimum at arbitrary high temperature 0. The canonical partition 
function demonstrates opposite behaviour: chiral symmetry is restored in all Z(N C ) 
phases and at the same temperature. We have shown this in our recent paper |13j by 
calculating the quark condensate in the CE in strong coupling SU (2) lattice gauge 
theory. 

It should be stressed that the canonical ensemble with respect to the total fermion 
number JTJJ also respects these features since it does not violate explicitly the Z(N C ) 



global symmetry. 

One may ask at this stage how to treat pure gauge theory where dynamical 
quarks are absent. The resolution is trivial if we remember that pure gauge theory 
is a limit of the full theory when the quark mass goes to infinity. This limit should 
be taken after the thermodynamical limit N —>■ oo, i.e. the pure gauge partition 
function is defined as follows 

1 N ° r 

Z PG = lim lim — y / T] dUiT] dVidVi e - Swi ~ ZkU ^- SF . (3) 

Thus, by definition a single quark does not appear in the spectrum and nonzero 
triality contributions are excluded by construction. Hence, the problem is only 
a careful mathematical treatment of the last equation. The puzzle mentioned in 
U and discussed in || does not appear. All the information can be extracted 
from correlators of PLs. In this respect the observation that the CE destroys the 
clusterization property has no sense as long as the PL itself has no meaning in 
the CE. If one wants to have this property and relation between the correlator in 
the limit R — > oo and the expectation value of the PL, the projection to another 
triality sector has to be performed. This later sector, however, does not present the 
ground state of QCD. < L > is a thermodynamical mixture of a heavy-light meson, 
a heavy- light baryon and so on (see for details, jnj). The correlator therefore gives 
the free energy of these states placed at distance R. 

Two questions are proper here. The first one concerns a suitable quantity which 
could characterize the critical behaviour of the system if one wants to use the CE. 
The second one is whether there is no contradiction between the screening of quarks 
and zero triality of the whole system. Indeed, it would naively seem that in the 
Debye screening phase there are no long-range forces. If we have a single quark at 
position r how can this quark feel what is going on in any other point far away 
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from r and, thus, how may the system have zero triality? It turns out that this 
later question is closely connected to the problem of detecting Z(N C ) charges at 
long distances. All these problems are the object of the next two sections. 



2 Order parameter for Z(N C ) charge and domain 
walls 

As an ensemble canonical with respect to N-ality does not break Z(iV c )-symmetry 
explicitly there appears the interesting question if Z(iV c )-symmetry can be broken 
spontaneously in this ensemble and which quantity could be used to detect nontrivial 
Z(N C ) charges at long distances. An operator which detects a nontrivial phase of 
some state in a region enclosed by a surface S is |L5| 



o - < mwn > m 

A[h > C) - <F(Z)xW(C)>> (4) 

and the limits E, C — > oo have to be taken. W(C) is a space-time Wilson loop 
nontrivial on the Z(N C ) subgroup. S may be thought as a set of plaquettes from 
one time-slice frustrated by Z(N C ) singular transformations in such a way that the 
dual links form a closed two-dimensional surface for this time slice. £ has to be 
chosen to enclose a line of the Wilson loop C going in time direction. The action 
of this order parameter is based on the Aharonov-Bohm effect: despite the absence 
of chromoelectric fields a singular potential influences particles at arbitrary long 
distances giving them a nontrivial phase during winding around the solenoid. In 
our case a Z(N C ) charge may be viewed as a kind of such singular solenoid placed 
at some space position. To get this configuration we should take the Wilson loop C 
and consider the limit C — > oo. Since the Wilson loop is a source of Z(N C ) charge we 
have the nontrivial charge and probe it via the "device" F(Y0). A nontrivial Z(N C ) 
charges may be detected only in the case when the fermionic screening is suppressed 
relatively to the Debye screening 



lim A(E, C) = exp 

S,C — >oo 



(5) 



where if (£, C) is the linking number of the surface £ and loop C. If the triality 
charge is totally screened one gets lim E) c_>oo A(£, C) = 1. 

Since the fundamental PL is a source of Z(N) charge at finite temperature we 
can insert the correlation function of PLs into @) instead of the Wilson loop 

A(T m - < mLpLn > (6) 

A ^ R) ~ <FM><LoL R > (6) 
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Another interpretation for the influence of F(E) is the following. F(E) tries to 
implement a phase change at the spatial surface E. If Z(N C ) symmetry is sponta- 
neously broken F(E) produces a stable interface between the volume enclosed by 
E and the surrounding vacuum. The different phases can be detected by a Wilson 
loop and the order parameter A(E, C). If the symmetry is unbroken -F(S) cannot 
produce a stable interface and the value of the Wilson loop is not influenced by 
.F(E) which is far apart. In the case of the grand canonical ensemble Z(N C ) symme- 
try is explicitly broken, there exists a preferred phase (zero) and therefore a phase 
change is annihilated after a short distance from E. At high and low temperatures 
lims,c^oo ^4(E, C) — 1 and cannot be used as an order parameter. Z(N C ) charges 
are completely screened by the nonzero iV c -ality contributions to the partition func- 
tion which are present in the grand canonical ensemble even at high temperature 
and in thermodynamical limit. As described above we expect however a nontrivial 
behaviour of this order parameter in the ensemble canonical with respect to 7V c -ality. 

To illustrate these ideas we studied a simple theory of Z(2) gauge spins coupled 
to massless naive fermions described in the CE by the partition function 

Z = \H H f d* x dt> x e s ™ +s *. (7) 

1 k=±ls t =±l J 

Sum over k is the sum over two iV c -ality sectors. We have denoted 

Sw = A ^ S p , (8) 

v 

where S p is a product of Z(2) gauge link variables s; around plaquette and 
1 1 

S F = a n E Sn(a:)[*x*x+n " + E S (x)k[V x V x+0 - ^ x+0 ^x}- 

" x,n x 

We fix a static gauge where all gauge spins in time direction s (x) are set to 1 except 
for one time-slice which includes the frustrated plaquettes E. We consider the region 
A > 1 and a<l which corresponds to the weak coupling high temperature phase. 
The fermionic screening is suppressed as a 2Nt . The Debye screening comes from the 
pure gluonic action and produces the formula for the correlation of PLs 

In < LqLr >oc exp[— MpR], (9) 

where Mr, is the Debye mass. It is straightforward to calculate < F(E) > in leading 
order in a 

< F(E) >= expf-^s + 0(a% (10) 
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where 5 « 2Q ^ff A . The main contribution comes from configurations of the gauge 
fields so(x) flipped in a volume enclosed by £ relatively to so(x) outside of this 
volume. Since the PL L in the origin penetrates this volume one time L changes 
its sign. The corresponding expansions in large A and in small a are converging. 
Therefore, in the limit of infinite E and R all the corrections go to zero and we find 

A(E, R) — —I. (11) 

In the region a>l the fermionic screening can dominate < LqLr > and we find 
A(E, R) = I. It follows that Z(2) may be broken in this region. We have thus shown 
that A(E, R) can be used in the CE as an order parameter to distinguish the phase 
of broken Z(N C ) symmetry from the phase where the triality charge can be detected 
at long range. The operator < F(E) > introduces a Z(N C ) string which becomes 
the boundary of a domain wall when the Z(N C ) symmetry is spontaneously broken. 
This domain wall is stable in pure gauge theory but appears to be unstable in full 
QCD if we treat it in the GCE. In the CE the domain wall is stable just as in pure 
gauge theory. Thus, the question of the realization of the domain walls in QCD is 
only the question whether the local Z(N C ) symmetry can be spontaneously broken 
in the CE in the sense described here. Application of these results to SU (N c ) theory 
is conceptually straightforward. Calculations however are much more complicated 
and could be advanced using Monte-Carlo simulations. 



3 Aq condensate in the CE 

We consider in this section another quantity which is a good candidate for character- 
izing the deconfinement phase, the so-called A condensate. It carries other useful 
features than the order parameter proposed above. Let us perform the following 
decomposition of the SU(N C ) PL V x before tracing (L x = TrV x ) 

V x = Z X V X: (12) 

where V x G SU(N C )/Z(N C ) and Z x e Z(N C ). For the invariant measure we have in 
this case 

w) = ^z/^(n (is) 

c z 

where D/j,(V) is the invariant measure on SU(N C )/Z(N C ) group. Let us recall, 
that the invariant measure on the SU(N C )/Z(N C ) group coincides with the SU(N C ) 
measure up to the restriction 

- - < arg[TrV] < -. (14) 
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We chose the static diagonal gauge for V x gauge field matrices. Then, V x has the 
diagonal representation with angles <pk obeying Eq. fll4l) . The corresponding con- 
tinuum potential can be written in the form /3g(A + a (x)) with a constant part A 
taking values as in ([14]) and a quantum fluctuating part a Q (x). This constant value 



may be interpreted as the A condensate, see |L7j for a review on this question. Our 
definition is, however, slightly different just because of the restriction ([14]). We have 
found that this definition is more proper both for the continuum theory and for the 
lattice theory. Z configurations are responsible for the disorder in the confinement 
phase and A takes values on the edge of the integration region (e.g., /3gA = ±ir 
for SU(2)). In the high temperature phase, it is possible to show that an actual 
value of Aq is shifted from the edge of the integration region forming a non-trivial 
saddle point in the SU(N C )/Z(N C ) subgroup [|IBJ (truly, we did not prove that this 
saddle point survives a transition to the continuum)^. 

We, thus, adjust the following definitions for the A condensate: 
1. On the lattice we define it as a saddle point configuration for the invariant 
integrals over the zero component gauge field matrix. Any such saddle point, if it 
exists, is invariant by the lattice construction. A physical meaning of this definition 
follows from the Hamiltonian formulation of lattice QCD: a constant saddle point is 
an (imaginary) chemical potential for the global colour charge |18[ . The existence of 
this constant saddle point configuration was demonstrated in [|T6j] for the Lagrangian 



formulation and in |T9[ for the Hamiltonian formulation. 

2. In the continuum the Aq condensate is defined as the position of the min- 
imum of the effective potential obtained by perturbative integration over gauge 
potentials over the vacuum Ai = 0, A =constant with fixed gauge. The proof of 
an independence of the condensation phenomenon of the gauge fixing as well as the 
independence of the physical quantities of the gauge in the theory with nonzero con- 
densate was done in Refs. |[20|| . Besides the same physical interpretation as before, 
the condensate plays the role of an infrared regulator of the theory. It was shown, for 
example, that there are no infrared divergences in the two-loop polarization tensor 



in a theory with such a condensate 2~I]. 

One important advantage of our definition of A Q is that it is insensitive to the 
Z(N C ) symmetry because V x is invariant under Z(N C ) transformation. Hence, it may 
be used both in a pure gauge theory and in a theory with dynamical fermions in both 



2 It should be mentioned that there is no a common opinion in the literature concerning the Aq 
condensate. Very different points of view have been stated and advocated during last five years, 
see fjj], [l6| and references therein. We are not going to discuss here all the issues on this subject 
conjecturing that a generation of Aq is possible both on the lattice and in the continuum space- 
time. We claim, however, that the value of Aq should be calculated after a summation over Z(N) 
configurations contained into the PL fill . 
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ensembles to reveal some features of the high-temperature phase. We think that Aq 
could be used as a sort of "order parameter" to distinguish the confining from the 
non-confining phase in the CE. Presumably, it is not an order parameter in the 
exact meaning, rather it shows preferred values of the SU(N C )/Z(N C ) configurations 
of the PL in different phases. In fact, this is very close to what the adjoint PL 
demonstrates in the 5*0(3) LGT at finite temperature [22]. Even though there is 



no global Z(2) symmetry in the 5*0(3) model, the phases of the adjoint PL behave 
differently in the confinement and in the deconfinement phases as described above. 
The 50(3) model was considered in this respect also in [|17|. We have proved for 
the simple example of a pure 50(3) model in the strong coupling region that the 
phases of the adjoint PL can indeed be used to test the deconfinement phase. In 
the continuum theory in the CE the A Q condensate was calculated in [^3| up to 



second loop order. It shows the expected behaviour at high temperatures. The loop 
expansion, however does not allow to observe a critical behaviour of the system. 
This is a question we address here. We are able now to demonstrate this expected 
behaviour in the strong coupling SU(2) LGT with Kogut-Susskind fermions. We 
work in the Hamiltonian formulation. The Hamiltonian of the lattice gluodynamics 
in the strong coupling approximation includes only the chromoelectric part 



^ links 



(15) 



where E{1) = id/d(A{) are the chromoelectric field operators. In this approach the 
chromomagnetic term can be treated perturbatively at g 2 — > oo. The determination 
of the partition function 

Z = Spexp(-f3H) (16) 

is connected with the summation over local gauge-invariant states. This is reflected 
by the symbol Sp in (|16"D, f3 is the inverse temperature. The corresponding physical 
Hilbert space is determined by Gauss' law. After the conventional procedure one 
gets the partition function of the form (see, for instance, flTi 



z= /n^)n 



(17) 



Ci is here the quadratic Casimir operator, 7 = 4^-. Vti is the character of Z-th 
irreducible representation of the SU{2) group 



sin(2/ + 1)0 



sin 1 



18) 
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Notice, that the invariant measure dn((f) x ) appeared after the representation of the 
Gauss' law delta-function on the 577(2) group. Following formulae (|T2l) -(|T4|) we 
represent SU{2) characters as 

= s 2 %{<f>). (19) 

Here, we may choose — | < < | which follows from ([T4]). We would like to 
calculate an effective potential for the phase which we set constant further. In 
this way, using a minimization procedure, we find a saddle point for integrals over 
the SU{2)/Z{2) subgroup. Substituting (|18|) and ([19]) into ( |TTD and summing up 
over all representations of the SU{2) group, one arrives to the effective potential of 
the form 

V eff {4>) = 2dlnF o (0) - (2d- l)lnsin 2 + ^l n Z /M (0), (20) 

where N s is a number of lattice sites and we have taken into account the measure 
contribution In sin 2 which is invariant under Z{2). Zim is the partition function 
of the Ising model 

ZlM= E n[ 1 + tanhyS - S -+n]Il[ 1 + McOS ^]' ( 21 ) 
s x =±l I x 



where 



and 



tanhy = £Sv < 22 > 



F o (0) = & 2 (e-\ 0) - Q 2 (e-\ 20), (23) 

F 1 (0) = e 3 (e^,O)-e 3 (e^,20). (24) 
Gj is the Jacobi theta-function. The quark contribution 

Z q = l[[l + M cos <f>8 x ] (25) 



was taken as in [21] with M = exp(— (3m). The first two terms from fl2Hp coincide 
exactly with the effective potential for pure the 5*0(3) model in the same approx- 



imation |]17 |. The existence of a nontrivial solution for can easily be established 
in this case for 7 > 7 C . For 7 < 7 C one finds = ±|. Thus, for this case our 



conclusion coincides with the one of Ref. [221. The third term in (EOT) establishes a 



difference from the 5*0(3) theory. It is clear, however, that this term cannot influ- 
ence drastically the qualitative picture in the deconfinement phase. We performed 
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numerical calculations to find maxima of the potential for the CE, where for the 
quark contribution we should take 

1 

Z T = E lit 1 + M cos (^ + k *)s x l (26) 

k=0 x 

whereas for the 50(3) model this contribution is 

z so(3) = J][i + Mcos 2 4 (27) 

x 

Detailed calculations with a chromomagnetic part of the Hamiltonian are at the 
moment in progress and will be presented elsewhere. The result for the effective 
potential ( pO|) including the quark contribution is, however, qualitatively the same 
both for the 50(3) model and for the SU(2) model in the CE, namely Aq forms 
a nontrivial saddle point above some critical value. Below this value <ft is if • 
Therefore, we established that A Q defined in such a manner, indeed can serve as 
an order parameter to probe phase transitions in a gauge theory with fundamental 
matter fields. 

As a last point we would like to mention an application of our result. We have 
already stressed that the Aq condensate has the meaning of an imaginary chemical 
potential for the global colour charge. Thus, if the condensate appears above some 
critical point, a spontaneous generation of nonzero colour charge might happen in 
the quark sector of the theory. This charge should have the same absolute value in 
all Z(N C ) phases. Since it must be zero in the confinement region this charge could 
be an exact order parameter for full QCD. If this phenomenon really happens it 
could lead to significant improvement of our knowledge of the quark-gluon plasma 
phase. 
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